STAR FORMATION IN THE NUCLEI OF SPIRAL GALAXIES 



T. Boker 1 2 

1 Space Telescope Science Institute - 3700 San Martin Drive - Baltimore, MD 21218 - U.S.A. 
2 Astrophysics Division - Space Science Department of ESA, ESTEC, Postbus 299, NL-2200 AG Noordwijk, The Netherlands 



ON 
ON 
ON 

o 

Q 

(N 



oo 

(N 

ON 
ON 

Oh 
i 

o 



Abstract 

Recent observations with the Hubble Space Telescope 
(HST) have revealed that a large fraction of late-type (Sc 
and later) spiral galaxies harbor a bright, compact stellar 
cluster in their dynamical centers. Statistics of the mass, 
age, and star formation history of these clusters as a func- 
tion of their host galaxy's Hubble type can be used to con- 
strain models of secular galaxy evolution. Since late-type 
spirals by definition do not possess a prominent bulge, 
their nuclear clusters are more easily separated from the 



The processes that lead to the formation of stellar clus- 
ters in the very center of galaxies are not well understood, 
and many questions remain to be answered: Is the cluster 
formation process self-regulating in the sense all clusters 
form over similar timescales? How stable are nuclear clus- 
ters? Is a stellar bar required to funnel gas into the central 
regions? 

The presence of such clusters is also interesting in the 
context of the ongoing debate on galaxy bulge formation. 



underl ying disk population. Their spectroscopic properties 



It is known from numerical simulations (e.g. Norman, Sell 
3d fc Hasan 199(f ) that the growth of a central mass 



can thus be studied from ground-based observations. Here, 
I will discuss plans for, and first results of, a program to 
study a sample of known nuclear clusters in late-type spi- 
rals. For one galaxy (IC 342), we have used high-resolution 
near infrared spectroscopy to determine the cluster mass 
directly via its stellar velocity dispersion. The analysis 
conclusively shows a very low mass-to-light ratio for the 
nuclear cluster in IC 342, indicative of a young cluster 
age (s=s 50 Myrs). From probability arguments, this result 
favors the scenario that such bursts are a recurrent phe- 
nomenon in late-type spiral nuclei. 
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1. Introduction 

Galactic nuclei are a unique physical environment as ev- 
idenced by massive black holes, active galactic nuclei of 
various flavors, and a star formation history that is very 
different from their host disks. While it is plausible that 
the nucleus might have unusual properties in galaxies with 
steep, cuspy light profiles, the gravitational potential van- 
ishes in diffuse, non-singular galaxies without a central 
black hole. In these cases it not at all intuitive why the 
nucleus should be the site of vigourous star formation. It 
therefore came as a surprise when recent HST data re- 
vealed the presence of prominent compact nuclear sources 



in many spiral galaxies of all Hubble types (Phillips et 



al. 1996: 


Matthews et al. 1999; 


Boker et al. 1999a; 


Car- 


olio 1999 


. More specifically, this is true for virtually all 



galaxies with exponential bulges - mostly galaxies of type 
Sc or later. Because the nuclear sources in most cases are 
resolved with HST, they are interpreted as stellar clusters. 



concentration can destroy stellar bars on short timescales. 
The bar undergoes a vertical "buckling" instability, dur- 
ing which stars are heated into orbits above the galaxy 
plane. The result is a triaxial system that much resembles 
a galaxy bulge. 

While this process is both plausible and appealing, 
there are also problems related to it. First, the observed 
high fraction of barred galaxies seems to contradict the no- 
tion that stellar bars are easily destroyed by central mass 
concentrations. Second, if the above scenario was true, 
then destruction of the bar would disable an effective fun- 
ncling mechanism for the gas to reach the galactic center. 
As a consequence, repeated nuclear cluster formation is 
prevented. 

Both these arguments can in principle be overcome 
if more and more massive bars form in a recurrent way. 
Building a galaxy bulge via repeated cycles of the bar for- 
mation - gas inflow - cluster formation - bar destruction 
cycle therefore seems like a possible alternative to theo- 
ries in which the bulge properties are largely determined 
during the collapse of the primordial gas cloud. 

Understanding the formation history of nuclear star 
clusters provides a new and important diagnostic to test 
the above scenario. Reliable statistics of nuclear cluster 
ages will directly test whether or not they form in a repet- 
itive cycle. For this reason, we have initiated a program 
to obtain high-quality spectra of the nuclear clusters of a 
sample of late-type, face-on spirals for which the cluster 
can be separated from the galaxy disk in ground-based 
observations. 

In this paper, I will describe our methods to derive the 
cluster ages from the observations and give an example for 
the successful application in Sec. 0. Sec. || discusses the 
preliminary evidence for repetitive cluster formation, and 
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Figure 1. Normalized CO-bandhead spectrum of the nu- 
clear cluster in IC 34-2, compared to four late-type tem- 
plate stars. All spectra are offset vertically for clarity; 
the abscissa is in fim. The heavy solid curve shows the 
best Gaussian-convolved template fit to the galaxy spec- 
trum, which has a = 33kms _1 . By contrast, the heavy 
dashed curve shows the ( unacceptable ) fit for a fixed value 
of a = 62kms~ 1 . The template spectra themselves show 
the instrumental dispersion of a — 5.5 km s -1 . 

Sec. |i| describes plans for an observational program which 
will put the analysis on much firmer statistical footing. 

2. Age-dating nuclar clusters 



It is w ell known from a number of population synthesis 
codes (e.g. Bruzual fc Chariot 1993| ; Lcithcrcr et al. 1999 ) 
that the mass-to-light ratio M/L of a stellar cluster is a 
sensitive age indicator. We have chosen the K-band for our 
analysis for the following reasons. Firstly, galactic nuclei 
are often affected by high and non-uniform dust obscura- 
tion. The near-infrared is much less sensitive to extinction 
than the optical: Ak ~ OA Ay - Secondly, the near- infrared 
emission is dominated by the evolved stellar population, 
and therefore gives a better representation of the stel- 
lar mass distribution than optical bands in which gaseous 
emission lines often contribute significantly. For both rea- 
sons, it is often easier to identify the true nucleus of a 
galaxy at near-infrared wavelengths. The third reason for 
choosing K-band is that it contains the 12 CO (0 — 2) band- 
head at 2.2936 fim, a prominent absorption feature which 
arises in the atmospheres of cool giants and supergiants. 
This feature is located in a spectral region with favorable 
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Figure 2. Relation between M/Lk (in solar units) and age 
(in years) for several model stellar clusters (from Lei- 
therer et al. 1999). Different curves are models that dif- 
fer in the slope a and upper-mass cut-off M up of their 
IMF. Solid: a = 2.35 and M up = 100 M Q ; long-dashed: 
a = 3.30 and M up = 100 M Q ; short-dashed: a — 2.35 
and M up — 30 Mq . All models assume solar metallicity 
(Z = 0.02), a lower-mass IMF cut-off at M low = 1 M , 
and an instantaneous burst. The dashed horizontal line in- 
dicates the value M/Lk = 0.047 inferred for the nuclear 
star cluster of IC 3^2. The models indicate that the cluster 
has log(age) < 7.7. 



observing conditions and has been shown to be a power- 
ful diagnostic for stellar k inematics in external galaxies, 
as discussed in detail by Gaffney, Lester fc Doppmanu 



(1995) 



Our methodology for deriving the M/L ratio is based 
on the "Gauss- Hcrmite Pi xel Fitting Software" developed 
by van der Marel (1994) . In brief, the observed surface 



brightness profile is deprojected under the assumption of 
spherical symmetry, and the three-dimensional velocity 
dispersion profile o~(r) is calculated by solving the Jeans 
equation for a spherical isotropic system. The results are 
then projected along the line of sight and convolved with 
the instrumental profile to yield a prediction for the ob- 
served velocity dispersion. This dispersion scales as a oc 
\J M/L, and we determine the M/L that produces the 
observed a. 

Two observational ingredients are needed for the above 
recipe to yield a reliable M/L value. In addition to a high- 
quality, high-resolution spectrum - yielding the cluster's 
velocity dispersion a - it is essential to also obtain an 
accurate surface brightness profile. Typical diameters of 
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Figure 3. Illustration of how both age (in years) and Ak of the nuclear cluster in IC 342 can be determined from 
the observed (V — K) b B = 4.57 and M/Lk — 0.047 x iQ- nA ( A K~OA5) _ Each panel corresponds to one of the three 
population models in Fig. [|. The dotted curves show the Ak for which the predicted V — K of the population model 
matches the observations, and the solid curves show the Ak for which the predicted M/Lk matches the observations. 
The intersections of the two curves indicate combinations of age and Ak that fit both observational constraints. The 
allowed combinations agree with the assumed extinction and the derived age, demonstrating self- consistency. 



nuclear clusters are 0.2" (Carollo 1999), therefore imaging 
with HST-like resolution is required. 

2.1. The showcase: IC 342 

We recently were able to successfull y apply the above 
method to the nuclear cluster of IC 342 ( Bokcr et al. 1999b| ). 



For these observations, we used the near-infrared spectro 
graph CSHELL ( |Tokunaga et al. 199C| ; [Green et al. 1993| ) 
at the NASA InfraRed Telescope Facility (IRTF) to de- 
termine the a from the shape of the CO (0-2) bandhead. 
Figure |l| shows the spectrum in comparison to a number 
of stellar templates that were used to derive a value of 
g = 33±3kms _1 . Together with an archival HST V-band 
image, we used the observed velocity dispersion to derive 
the cluster mass-to-light ratio according to the method 
outlined in H. Our data imply M/L= 0.05 Mq/L q and a 



2.2. Dealing with extinction 

In galactic nuclei, dust extinction can significantly affect 
the observations, even at nearinfrared wavelengths. How- 
ever, if values for both color and mass-to-light ratio are 
available from the observations, it is possible to perform 
a consistency check on the assumed extinction. This is 
demonstrated in Fig. || which for three different IMF's 
shows possible combinations of age and extinction which 
explain the observed (V — K) color and M/Lk ratio, re- 
spectively. Locations where the two lines intersect indi- 
cate age/extinction combination which are consistent with 
both observables. The allowed combinations have ages in 
the range 10 6 - 8-7 - 8 years, and Ak between 0.37 and 0.57. 
This range includes the value Ak ~ 0.45 (i.e., Ay = 4.0) 
that was adopted for our analysis based on the results of 
BFG97, and the age that we derived for the nuclear clus- 
ter. 



In order to derive the cluster age, we compared these 
results to predictions from population synthesis models 



such as the "Starburst99" package of Leitherer et al. (1999) 



We used an instantaneous burst model for our compari- 
son which is justified by the high equivalent width (EW) 
of the CO absorption as demonstrated in an earlier pa- 
per (Boker, Forster-Schreiber & Gcnzcl 1997, hereafter 
BFG97). Fig. ^| shows that - quite independent of the as- 
sumed Initial Mass Function (IMF) - the nuclear cluster 
in IC 342 is rather young, namely less than 50 Myrs. 



3. Repetitive cluster formation? 

Given the young age and low mass-to-light ratio of the 
IC 342 nuclear cluster, one is led to ask whether its for- 
mation was the first of this kind. The presence of an 
underlying older population is not ruled out by our ob- 
servations, but there are constraints on the amount of 
light and mass that such a population could contribute. 
These constraints are based on the observed EW of the 
CO bandhead at 2.29 ^m which is given by EW to tai = 
EW young - /(EW youl 



EW id). Here, / is the fraction 



of the total light contributed by an old population. The 
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EW for single-burst population models decreases with age, 
so the observed EW sets an upper limit on /, namely 
/ <0.14. 

Although an old population can only contribute a small 
fraction of the observed light, its contribution to the to- 
tal mass of the cluster may be significant. For the same 
mass, a burst that happened 10 9 5 yrs ago would be 3 
magnitudes fainter in the K-ba nd than the young popula - 
tion that dominates the light ( Bruzual fc Chariot 1993| ). 
An old population that formed 10 10 2 yrs ago, roughly the 
age of the universe, would be 5 magnitudes fainter. In the 
latter case it could be 15 times as massive as the observed 
young population and still only have / « 0.14. The total 
mass of the cluster is fixed, so in this scenario the young- 
population would need to have a mass-to-light ratio that 
is 15 times smaller than previously inferred, which implies 
M/Lk ~ 0.003. This is roughly the smallest value that 
can be plausibly explained with a single-burst population 
model (cf. Fig. ||), and cannot be ruled out. So while we 
can conclude that the light of the cluster is dominated 
by a young population, this is not necessarily true for the 
mass. 
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Figure 4-. Histogram of nuclear stellar cluster luminosities 
for the galaxy sample of Carollo (1999). Only galaxies with 
detected nuclear clusters are included. IC 342 is indicated 
by the arrow, but is itself not included in the sample used 
to construct the histogram. The labels at the top indicate 
log(age) for each bin, using a simple model that assumes 
that all clusters are identical except for their age. The lu- 
minosities were calculated from the stellar population mod- 
els of Bruzual & Chariot (1993), and are normalized to 
log(age) = 7.4 for IC 342. 



For the Galactic center, spectroscopy of individual stars 
is possible, and such data provide solid evidence for mul- 
tiple periods of star formation (Krabbe et al. 1995). In 



general, however, it is a difficult task to identify multiple 
star formation events from the integrated light of stellar 
populations, and even more so to find a unique solution 
to their star formation history. The best hope for answer- 
ing whether or not stellar clusters form repeatedly in the 
nuclei of spiral galaxies stems from a statistical analysis 
of their frequency and age. 

3.1. Statistics so far 



IC 342 is not the only galaxy known to harbor a young 
population of stars in its nucleus. The Milky Way, for 
example, underwent a recent starburst some 3 Myrs ago 
(|Krabbc et al. 1995|), and both M31 and M33 hav e blue 



nuclei that are quite possibly young star clusters (Lauei 
et al. 1998). If nuclear cluster formation is a Poissonian 
process with equal probability per unit time, one would 
expect to find clusters with ages in the range of 10 7±0 5 , 
l 8±o.5 j and 1Q 9±o.5 yearg in the ratio ! ■ io : 100. 

In order to test this prediction in the absence of reliable 
age estimates for a large number of clusters, we took the 
following approach. Under the very simplifying assump- 
tion that all nuclear clusters have the same mass as the 
IC 342 cluster, and formed in an instantaneous burst, one 
can estimate their ages directly from lumi nosities. This is 
shown in Fig. |] for all clusters found in the Carollo (1999) 
sample. While the interpretation of the Fig. |4| is certainly 
hindered by the extremely simplifying assumption of equal 
mass and formation history for all clusters, it nevertheless 
seems that there is an overabundance of relatively young 
clusters which can be interpreted as preliminary evidence 
for a repetitive process. 

4. Observational plans 

In order to improve the available dataset on nuclear clus- 
ters, we are planning the following observations: 

— High-resolution CO bandhead spectroscopy (NIRSPEC 
@ Keck - approved) of all galaxies found to date with 
HST to harbor a bright (K < 16) nuclear cluster with 
the goal to derive stellar velocity dispersions. 

— Optical Echelle spectroscopy (BC spectrograph @ Stew- 
ard 90" - data taken) of the same targets to constrain 
their stellar populations. Figure || shows "quick- look" 
results of the data reduction for two of the sample 
galaxies. In both cases, the (disk-subtracted) cluster 
spectra show deep Balmer absorption lines indicative 
of fairly young cluster ages. The detailed analysis of 
these spectra is ongoing, and will be published sepa- 



rately (Mazzuca et al. 2000) 



NIR imaging of about 100 late-type, face-on spirals 
(MAGIC @ MPIA 2.2 m- approved) in order to iden- 
tify a larger sample of prominent nuclear cluster. 
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Figure 5. Echelle spectra of the nuclear clusters of NGC 4449 (top) and NGC 4559 (bottom). The spectrum of the 
underlying disk/bulge population has been subtracted. The feature at 5577 A in the spectrum of NGC 4559 is a residual 
skyline. Both clusters show deep Balmer absorption features, indicative of a young cluster age. 



- A cycle 9 HST WFPC2 R-band snapshot survey (pro- 
posed) to identify a large number of nuclear clusters 
for follow-up spectroscopy and obtain their surface 
brightness profiles. 
The combination of ground-based spectroscopy and HST 
imaging will allow us to measure cluster masses and ages in 
a way similar to that described for IC 342 for a large sam- 
ple of galaxies. The cluster mass and age distribution as a 
function of host galaxy Hubble type will answer whether 
the formation of the nuclear clusters in spiral galaxies is 
a recurring phenomenon or not, and how - if at all - these 
clusters are related to the secular evolution of the expo- 
nential bulges in which they reside. 
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